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Abstract 

The hyperfine splittings in heavy quarkonia are studied using new 
experimental data on the di-electron widths. The smearing of the spin- 
spin interaction is taken into account, while the radius of smearing is 
fixed by the known J/ijj — Vc{^S) '0(2»S') — v'd'^S) splittings and 
appears to be small, Vss ~ 0.06 fm. Nevertheless, even with such a 
small radius an essential suppression of the hyperfine splittings (~ 
50%) is observed in bottomonium. For the nS bb states (n = 1, 2, ...6) 
the values we predict (in MeV) are 28, 12, 10, 6, 6, and 3, respectively. 
In single-channel approximation for the 3S and 45 charmonium states 
the splittings 16(2) MeV and 12(4) MeV are obtained. 

1 Introduction 

At present two spin-singlet S-wave states rjc^lS) and ric{2S) are discovered 
[lJ-[3 ; still, no spin-singlet rih{nS) levels have been seen jl]. Though the 
masses of the rih{nS) were predicted in many papers the calculated 

hyperfine (HF) splittings, Ahf('^'S') = M{n^Si) — M{n^So), vary in a wide 
range: from 35 MeV up to 100 MeV for the bb IS state and for the 25* 
state between 19 MeV and 44 MeV |llj. However, at the modern level of 
the theory and experiment there exist well-established limits on the factors 
which determine the spin-spin potential Vhf(^) in heavy quarkonia. First of 
all, the wave function (w.f.) at the origin for a given n ^Si{cc or bb) state 
can be extracted from di-electron width which are now measured with high 
accuracy jT2l CB]- Concerning the quark masses, the pole (current) mass. 



present in a correct relativistic approach, and the constituent mass, used in 
nonrelativistic or in more refined approximations, are also known with good 
accuracy ^3 E] • Therefore the only uncertainties comes from two sources: 

First, in perturbative QCD there is no strict prescription how to choose 
the renormalization scale in the strong coupling anp, entering V^F^r). 

Secondly, the role of smearing of the spin-spin interaction is not fully 
understood and the true size of the smearing radius r^g is still not fixed. 

Moreover, the masses of higher triplet and singlet states can be strongly 
affected by open channel(s), thus modifying the HF splittings. 

In our calculations the smearing radius Tss is taken to fit the J /ip — ric{lS) 
and ipi'^S) — ric{2S) splittings. To reach agreement with experiment it is 
shown to be small, r^s < 0.06 fm. Our value r^^ = 0.057 fm practically coin- 
cides with the number used in Ref . ^U] . However, in spite of this coincidence 
the splitting Ai = T(15') — rjb{lS) = 28 MeV in our calculations appears to 
be two times smaller than that in Ref. JUl; where Ai = 60 MeV. 

From our point of view the use of the w.f. at the origin \Rn{0)\'^^p, ex- 
tracted from di-electron widths, is the most promising one, because these 
w.f. take implicitly into account the relativistic corrections as well as the 
infiuence of open channel(s), in this way drastically simplifying the theoret- 
ical analysis. A comparison of these w.f. with those calculated in different 
models puts serious restrictions on the static potential used and also on 
many-channel models. 

The HF interaction is considered here in two cases. First one corresponds 
to the standard perturbative (P) spin-spin interaction with a 5-function: 

VFsir) = si ■ S2^as{fi) (l + 5(r) = Si ■ S2VnF{r), (1) 



which in one-loop approximation gives the following HF splitting |6j: 

(2) 



where p = ^/3o — | — | ln2 and the second term in brackets is small: < 0.5% 
in bottomonium [nj = 5) and < 3% in charmonium [nj = 4). 

It is very probable that 6{r) may be considered as a limiting case and the 
"physical" spin-spin interaction is smeared with a still unknown "smearing" 
radius. For the Gaussian smearing function 

(5(r) ^ ^ / r2rfrexp(-/5V2) (3) 
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the splitting can be rewritten as 

^UnS) = l^UmMO)\' (l + -p) , (4) 
9 uJq \ it J 

where by definition the "smearing factor" C,n{(3) is 

UP) = ■^J^^2 J \Rn{r)\'exp{-P'r'ydr. (5) 

The general expression is evidently kept for any other smearing prescrip- 
tion which may differ from Eq. (jSj). 

2 Wave function at the origin 

The w.f. at the origin is very sensitive to the form and parameters of the 
gluon-exchange interaction and also to the value of the quark mass used. 
Therefore we make the following remarks: 

a. To minimize the uncertainties in the w.f. at the origin, Rn{0), we shall 

use the w.f. extracted from the experimental data on leptonic widths 
and denote them as |-Rn(0)|exp- In this way the relativistic corrections 
to the w.f. and the infiuence of open channel(s) are implicitly taken 
into account. 

b. In Eqs. (0) and (jH) the constituent mass Ug enters (this fact can be 
rigorously deduced from relativistic calculations O CH]): ^^q{nS) = 
{\f^~+^f^)ns, where under the square-root the pole mass mq = mq 
(pole) is present. This mass is known with good accuracy and we take 
here mb(pole) = 4.8 ±0.1 GeV and mc(pole) = 1.42 ±0.03 GeV, which 
correspond to the well-established current masses mi,{rhf,) = 4.3(1) 
GeV, fhcijhc) = 1-2(1) GeV j3j, while the constituent masses he also in 
a rather narrow range for all nS states, both in charmonium and bot- 
tomonium: UbinS) = 5.05 ± 0.15 GeV, LudnS) = 1.71 ± 0.03 GeV(n > 
2), and udlS) = 1.62 ± 0.04 GeV(n = 1). Note, that just these mass 
values are mostly used in nonrelativistic calculations, thus implicitly 
taking into account relativistic corrections. 



3 



c. The leptonic width of the n ^Si states in heavy quarkonia are defined 
by the Van Royen-Weisskopf formula with QCD correction 7q, 

Tee{n 35i)|exp = |^n(0) |Lp7Q, (6) 

n 

where eg = | for a 6(c) quark, a = (137)^\ M„ = M{n ^Si), and 
•jQ^nS) = 1 — |^as(2mQ), with the renormahzation scale in ag equal to 
2mQ (pole), as in Refs. 01101 ^^^^ Vb 77 decay pTOI- In some cases 
/i = Mn is also taken, but with accuracy < 1% both choices coincide (here 
2mh = 9.6 GeV and 2mc = 2.9 GeV are taken). 

Since for = 5 the QCD constant is well known from high energy 
experiments the factor •jf, is also defined with a good accuracy. For 
A^{3-loop) = 210(10) MeV, which corresponds to aj^{Mz) = 0.1185, one 
has 76 = 7b„ = 0.700(5) and as(2mb) = 0.177(3). In charmonium (n/ = 4) 
for = 0.260(10) MeV the coupling a,(2mc = 2.9 GeV) = 0.237(5), 



7c = 0.60(2). Then the w.f. at the origin, extracted from the di-electron 
width ©, 

|^„(0)|L„^^%i^, (7) 



implicitly takes into account the relativistic corrections as well as the infiu- 
ence of open channels, which gives rise to smaller values for |i?„(0)| as well 
as for the HF splitting. The extracted values of |-Rn(0)|gxp in the bb and cc 
systems are presented in Table [H 

The extracted |-Rn(0)|exp can be compared to the predicted values, which 
chiefiy depend on the strong coupling used in the gluon-exchange term. In 
particular, if the asymptotic-freedom behavior of astatic('^) is neglected, then 
theoretical numbers can be 2 — 1.5 times larger than |-Rn(0)|exp; even for the 
T{nS) (n = 1, 2, 3) states, which lie far below the BB threshold JS]. 

Here, as well as in our analysis of the spectra and fine structure splittings 
in heavy quarkonia [7l[T31[Tnj, we use the static potential Vsir) in which the 
strong couphng a_B(r) is defined as in background perturbation theory: 



iasir) 8 /• singr 

VB{r) = ar - , aB{r) = — dq 



Pi In tB 

Ml 



3 r PoJq tij(g) 



(8) 



where tsiq) = In ^2 ■ Here Mb = 0.95(5) GeV is the background mass, 
AB^nf) is expressed through Aj^{nf) and in 2-loop approximation KBirif = 
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Table 1: The w.f. |^„(0)|^^p (in GeV^) and the leptonic widths Tee{T{nS)) 
and Teeii^inS)) (in keV)"'^) (7^ = 0.70, 7^ = 0.60). 





bb 


cc 


Tee (^*S')exp 


-Rn(O) 


Tee (^*S')cxp 


-Rn(O) cxp 


1 c 
10 


1.314(29)") 


7.094(16) 

7010/1 p:\ 
( .ZioyiO) 


5.40(22) 

C ftO (OA \ 

0.0o(^z4 ) 


0.911(37) 
u.yoy (^4u 


2S 


0.576(24)") 
0.616(19)^) 


3.49(15) 
3.73(12) 


2.12(12) 
2.54(14) 


0.51(3) 
0.61(3) 


3S 


0.413(10)^*) 


2.67(7) 


0.75(1) 
0.89(8) 


0.22(1) 
0.26(2) 


AS 


0.25(3)") 


1.69(20) 


0.47(15) 
0.71(10) 


0.16(5) 
0.24(4) 


5S 


0.31(7)") 


2.21(49) 






6S 


0.13(3)") 


0.95(22) 







") The upper values of the leptonic widths in bottomonium are taken from 
PDG 13] and the lower values of ree(^5') are taken from the CLEO data fr2\ . 
'') The upper entries in charmonium are taken from [3] and the lower ones 
from [TH] . 



5 



Table 2: The factor Sn (9) for the potential VB{r) (8) in charmonium and 
bottomonium. 





1^ 


2S 


3^ 


4^ 


5^ 


6S 


bb 


1.08(4) 


1.02(4) 


1.02(4) 


0.72(9) 


1.03(22) 


0.47(10) 


cc 


1.01(4) 


0.82(5) 


0.41(2) 


0.32(10) 







4) = 360(10) MeV and Aeirif = 5) = 335(5) MeV HU; the string tension 
a = 0.18 GeV^. Our calculations show that in bottomonium (in single- 
channel approximation) the potential Vsir) gives values of |-Rn(0)| theory '^^ry 
close to the values |-Rn(0)|exp- For illustration in Table |21 the ratios 

„ _ I -Rn (0)1 cxp /f.\ 
l-'^K^J I theory 

are given for all known nS levels in charmonium and bottomonium. 

As seen from Tabled using potential VB{r) the influence of open channels 
in bottomonium appears to be important only for the 45" and 65* levels, while 
for the other states single-channel calculations are in good agreement with 
experiment. This is not so for many other potentials [S] and it means that 
any conclusions about the role of open channels cannot be separated from 
the QQ interaction used in a given theoretical approach. 

In charmonium the effect from open channels is much stronger and already 
reaches ~ 60% for the 35' and the 45* states {Sn ~ 0.4) and about 20% for 
the ip{2S) meson. 



3 Hyperfine splitting 

Now we discuss the HF splitting for both bottomonium and charmonium. 
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3.1 Bottomonium 



In bottomonium the HF splittings are considered in two cases: 

A. No smearing effect, i.e., in Eq. (0)) the smearing parameter ^ft„ = 1.0 (Vn). 

B. The smearing parameter (0) is calculated with (3 = \/T2 GeV, or a 

smearing radius Vgs = = 0.057 fm. 

Unfortunately, at present there is no a precise prescription how to choose 
the renormalization scale in the HF splitting (2): in aj^{jl) the scale /i = 
m6(pole) ^ 4.80 ± 0.01 GeV is often used. With Ajrslrif = 5) = 210(10) 
MeV (just the same as in our calculations of % (6)) one finds 

a,{bb, p) = aj^{4:.8 GeV) = 0.21(1). (10) 

With this a,(/i) and {RnWe^^ from Table HJ one obtains the HF splittings in 
bottomonium presented in Table El second column. (The numbers in Table 
3 contain experimental errors coming from ree(^>S') P] (first number) and 
theoretical errors (second number)). For the T{nS) states {n = 1,2,3) the 
calculated HF splittings (^„ = 1.0) appear to be very close to the splittings 
from Refs. 0. 

If smearing of the HF interaction (3) is taken into account (the smearing 
radius, r^s = P'^ = 0.057 fm for /? = a/12 GeV, is taken to fit the experi- 
mental values of the J / ip — rjcilS) and ip{2S) — rjci^S) splittings), then even 
for such a small radius Ahf('T"S') turn out to be 50% {n = 1, ... ,4), 60% 
{n = 5, 6) smaller as compared to the "nonsmearing" case. In particular, the 
T(lS)-r/b(15) splitting turns out to be 28 MeV instead of 51(4) MeV for 

= 1.0. For higher excitations very small splittings, Ahf ~ 6 MeV and 3 
MeV for the 55* and 65 states, are obtained, see Table El 

Note that our value of rss = 0.057 fm is very close to that from Ref. [Tllj 
where Tss = 0.060 fm is taken. However, in spite of this coincidence our 
numbers are about two times smaller than in ^], where T(IS') — r]b{lS) 
= 60 MeV is obtained. For the 2S state our value of the splitting is 12 
MeV, still smaller than 20 MeV in pUj. From this analysis it is clear that 
the observation of an rib{nS) meson could clarify the role of smearing in the 
spin-spin interaction between a heavy quark and antiquark. 
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Table 3: Ayp(nS') (in MeV) in bottomonium for <yj^{jj,) = 0.21, cjf, = 5.10 
GeV and |^n(0)P from Tabled 





6 = 1-0 
(no smearing) 


^bn for the smeared HF interaction 
with (3 = ^12 GeV, 
rss = 0.057 fm 


IS 


51(4) (4) 


28(2) (3) 


2S 


25(3) (2) 


12(2) (1) 


3S 


22(5)(2) 


10(2) (1) 


4S 


12(3) (1) 


5.1(2) (1) 


5S 


16(2) (1) 


6.4(1) (1) 


6S 


7(2) (1) 


2.7(1) (1) 
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3.2 Charmonium 



Also in charmonium the sphtting (jlj) in fact depends on the product as{fi)-C,n, 
therefore it is convenient to discuss an effective HF couphng: aupinS) = 
ots{jln)icn^ which is the only unknown factor. (The masses udnS) may be 
specified for different nS states [14,15].) 

As discussed in j7j, the experimental splittings J /ip — rjc (IS) and ifj (2S) 
— '/7c(2S) can be fitted if different values of anr for the IS and 2S states are 
taken namely, aHF(lS) ~ 0.36 and q;hf(2S) ~ 0.30. Such a choice implies 
two possibilities. The first one, case A, is 

A. as(Aii) = 0.36, as(/U2) = 0.30, ^^(/is) = as(/i4) < 0.30, icn = 1-0 (Vn), 

(11) 

i.e., the renormalization scale is supposed to grow for larger excitations. In 
particular, for A^^(2 — loop) = 270 MeV one finds /ii = 1.25 GeV ~ fhc{fhc) 
while the scale ^2 = 1-60 GeV is essentially larger. For this choice of anr 
the perturbative HF splittings are given in Table 01 second column. 

Besides, we have also calculated the contributions coming from the NP 
spin-spin interaction. In bottomonium their values are small, A^p(nS') < 1 
MeV, and can be neglected. In charmonium, as well as in light mesons, the 
situation is different, e.g. due to the NP spin-spin interaction in the IP cc 
state a cancellation of perturbative and NP terms takes place P^HJ. As a 
result, the mass difference Mcog(xcj) — M{hc) = (+1 ± 1) MeV turns out to 
be close to zero or even positive, in accord with experiment PH]- The values 
of A2p(n5') are given in Table IH fourth column. 

Thus one can conclude that in case A with different renormalization scales 
fin, the splittings J/ip — Vci^S) ^^'^ 'V^(25') — rjc{2S) can be obtained easily 
in agreement with experiment. 

If the renormalization scales supposed to be (almost) equal for all 

nS states: 

B. a,(/i„ ^ m, = 1.25 GeV) = 0.36, (12) 

then to explain the relatively small ilj{2S)—ric{2S) splitting, a smearing effect 
needs to be introduced. Then for the potential used, the values ^„(cc) = 
0.84, 0.80, 0.78, and 0.76 for the 15", 25", 35", and 45* states, respectively, 
are calculated. In this case the A^p{cc,nS) are also given in Table lU For 
the higher 35" (45") levels our predicted numbers are about 21(15) MeV (no 
smearing) and 16(12) MeV (with smearing), i.e., the difference between cases 
A and B is only ~ 20%. Notice that in case B the NP contribution improves 
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Table 4: The splittings A^p(?T,S') and A2p(r;,S') (in MeV) in charmonium"^ 





(no smearing: C,c = 1-0) 
asifii) = 0.36;as{fJ^n) = 0.30 
(n = 2,3,4) 


r,, = 0.29 GeV-i 

«s(yUn) = 0.36 

(n = l,...,4) 


ANP(n5)'') 
G2 = 0.043 
GeV^ 


1^ 


117(5) 


102(6)'^) 
108(7)") 


9±2 


experiment 

T Inly r, H 


117(2) 


117(2) 




9 Q 




55(4) 


0.0 It i.O 


experiment 
^l; (2S) -r]c (2S) 


48(4) 


48(4) 




35 


21(2) 


16(2) 


2± 1 


4S 


15(4) 


12(4) 


1.5 ±0.5 



The w.f. |-R„(0)p, taken from Table 1, correspond to ree(^>S') from PDG[3]. 
^) The NP splittings are calculated in [17]. 

Here |^i(0)|2^p = 0.959 GeV^ and |^2(0)P = O.eiGeV^ from the CLEO 
data [13] are taken. 
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the agreement with experiment for J/ijj — rjcilS). As a whole, in charmonium 
the smearing effect appears to be less prominent than in bottomonium. 

4 Conclusions 

Thus we come to the following conclusions: 

1. In bottomonium AHp(n5') appears to be very sensitive to the smearing 
of the spin-spin interaction. Due to this effect the sphtting decreases 
from 51 MeV to 28 MeV for the IS" state and from 25 MeV to 12 MeV 
for the 25* state; very small values are obtained for higher states. 

2. In charmonium there are two possibihties to describe Ahf(I'S') and 
Ahf(25'), which are known from experiment. First one refers to a 
different choice of the renormalization scale: /ii =1.25 GeV and ~ 
1.60 GeV for the l-S" and 25" states, if the smearing effect is absent. The 
second possibility implies equal renormalization scales /i„(n = 1, . . . , 4) 
for all nS states. Then to explain the ip{2S) — ric{2S) sphtting the 
smearing of the spin-spin interaction needs be taken into account. We 
also expect that for the IS level a small contribution (~ 8 MeV) comes 
from the NP spin-spin interaction. 

3. The ipi'iS) — rjci'iS) splitting in single-channel approximation is pre- 
dicted to be around 16(2) MeV, without and 12(4) MeV with smearing 
effect. 

In order to understand the true role of the smearing effect in the spin-spin 
interaction the observation of an r]b{nS) is crucially important. 
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